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Abstract 

It is shown that for a relativistic particle moving in an electromagnetic 
field its equations of motion written in a form of the second law of Newton 
can be reduced with the help of elementary operations to the Hamilton- 
Jacobi equation. The derivation is based on a possibility of transforming 
the equation of motion to a completely antisymmetric form. Moreover, by 
perturbing the Hamilton-Jacobi equation we obtain the principle of least 
action. 

The analogous procedure is easily extended to a general relativistic 
motion of a charged relativistic particle in an electromagnetic field. It 
sis also shown that the special-relativistic Hamilton-Jacobi equation for 
a free particle allows one to easily demonstrate the wave-particle duality 
inherent to this equation and, in addition, to obtain the operators of 
the four-momentum whose eigenvalues are the classical four-momentum 
03.20. +i, 03.30. +p 

In analytical mechanics we arrive at Newton's second law (the experimentally 
verified phenomenological equations with the observable parameters) by postu- 
lating the principle of least action. In turn, the action S obeys the Hamilton- 
Jacobi equation. The latter is a partial differential equation of the first order. 
A transition from Newton's second law to the Hamilton-Jacobi equation can be 
achieved with the help of the algorithm for transforming a system of ordinary 
differential equations into a partial differential equation. Despite the fact that 
such transformation algorithm is well-known (e.g., pQ) the actual transformation 
of the equations of motion of a charged relativistic particle in the electromag- 
netic field into a respective PDE (the Hamilton-Jacobi equation) is not quoted 
in the physical literature to the best of our knowledge. The usual approach to 
the problem of derivation of the Hamilton-Jacobi equation is to heuristically 
introduce classical action S and to vary it (for fixed initial and final times). 



'Department of Physics, University of the Pacific, Stockton, CA. 95211; E- 
mail: agranik@uop.edu 



1 



Here we provide an elementary derivation of the Hamilton-Jacobi where the 
concept of action emerges in a natural way by considering the momentum as 
a function of both temporal and spatial coordinates. This can be seen by con- 
sidering first a non-relativistic classical particle moving from p. A to p.F (see 
Fig.l). The particle can do that by taking any possible paths connecting these 
two points. Therefore for any fixed moment of time, say t = 1 the momentum 
would depend on the spatial coordinate, that is p = p{x, t). In a sense we have 
replaced watching the particle evolution in time by watching the evolution of 
its velocity (momentum) in space and time. This situation is analogous to the 
Euler's description of motion of a fluid (an alternative to the Lagrange descrip- 
tion). The other way to look at that is to consider a "flow" of an "elemental" 
path and describe its " motion" in terms of its coordinates and velocity (deter- 
mined by a slope of the path at a given point). 

This allows us to represent Newton's second law for a particle moving in a 
conservative field U (x) as follows 

|=f+Vv)p = -vc/ (1) 

at at m 
We apply curl to both sides of this equation and get 

dp d 1 -* 

curl — = — curlp -\ curllp ■ V)p = (2) 

at at m 

Using the vector identity 

- Va 2 

[a ■ V)a = — h curia x a (3) 

we rewrite Eq.JSJ) 

—curlp H curUcurlp x p) = (4) 

at m 

One obvious solution to Eq.Q is 

curlp = 

similar to an irrotational motion in Euler's picture of a fluid motion. 
Eq.QJ implies 

p = VS (5) 

where S(x, t) is some scalar function. Generally speaking, we can choose the 
negative value of VS. The conventional choice is connected with the fact that 
the corresponding value of the kinetic energy has to be positive. Upon substi- 
tution of Eq.0 into Eq.0J we obtain with the help of Eq.Q the following 
equation 

v <§ + i< vs > 2 + c, > = » < 6 > 
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Figure 1 : A few paths of a path set connecting the initial and the final points travelled 
by a particle in t = 3sec. ft is clearly seen that particle's velocity ( momentum) is a 
function of both coordinate x and time t 



In turn Eq.© means that 

^ + -(vs)» + cr = /( t ) (7) 

where f(t) is some function of time. By introducing a new function 

S' = S- J f{t)dt 

we get from (JBJ) the Hamilton- Jacobi equation with respect to the function S 1 
(representing the classical action): 

^ + _L (VS <)* + t / = (8) 



If we use relation JSJ, v — p/m, and drop the prime at S' , the Hamilton- Jacobi 
equation can be rewritten as follows 

OS _ mv 2 , , 

- + e.vs = — -u (9) 

Since by definition 

dS dS _ „_ 

we obtain from © the expression for the action S by integrating |([JJ from p. A 
to p.F 

f tp mv 2 f tp 
S= (— U)dt = j L{S,v,t)dt (10) 

Jt A 2 ^tA 
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where L(x, v, t) — mv 2 /2 — U is the lagrangian of a particle of mass m. 

Now we can arrive at the principle of least action (without postulating it a 
priori) directly from the Hamilton-Jacobi equation. To this end we subject the 
action S to small perturbations SS <C S and (by dropping the term (\7SS) 2 ) 
get from JHJ the equation with respect to SS 

^ + i(VS).(V^)=0 (11) 
Since 'VS/m = v Ea. (|ll(l represents the substantial derivative of SS: 

This means that 

SS = const (13) 

Thus on one hand, for a specific function S satisfying the Hamilton-Jacobi equa- 
tion the respective perturbations SS = const. On the other hand, according to 
Eq. (|10fl the action S is defined on a set of all possible paths connecting point 
A and point F. This means that perturbations SS correspond to perturbations 
of all these path. 

However, only for one of these paths SS — const, according to l|13|l . To deter- 
mine this constant we take into account that at the fixed points A and F the 
paths are also fixed, that is the respective perturbations SS — at these points. 
Therefore only for the specific path determined by the Hamilton-Jacobi equa- 
tion ( that is by the second law of Newton) SS — 0, thus yielding the principle 
of least action: 

S I L(x,v,t) = (14) 
Jt A 



The above derivation serves as a guide for a derivation of the Hamilton-Jacobi 
equation for a relativistic particle of charge q and mass m moving in the electro- 
magnetic field. Our approach is to reduce the respective equations of motion to 
the form which would be analogous to an irrotational motion in Euler's picture. 
The very structure of the space-time metric allows one to arrive at the required 
result in a natural way. 

Therefore we begin with the second law of Newton for a relativistic charged 
particle of a charge q and mass m moving in the electromagnetic field: 

^ = q[E a + e a ^v^B^} (15) 

where Greek indices a, (3, 7, ... take the values 1, 2, 3, e a ^ J is the absolutely an- 
tisymmetric tensor of the third rank, p a = mv a /(l — v s v s ) 1 / 2 is the momentum 
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of the particle, E a is the electric field , v a = v is the velocity of the particle and 
B a is the magnetic field. 

For the subsequent analysis we cast Eq. i|15fl into the standard co - and contra- 
variant forms. To this end we use the metric g tk = gik = [1, — 1, — 1, —1] and 
use units where the speed of light is c = 1. In this metric x° = xq = t, 
x a — x — ~ x a , the four- potential A 1 (A , A a ) whose scalar part A = <fi (where 
4> is the scalar potential) and A a = A is the vector potential, and the roman 
indices i,j,k ,. . . take the values 0,1,2,3,. From the Maxwell equations then 
follows (e.g.0) that the electric field E a intensity and the magnetic induction 
B a are 

E ~- { dx^ + ^ ] (16) 

Ba = ^ J d^ (17) 

Using (|17(l we express the second term on the right-hand side of Eq. I|15|) in 
terms of the vector-potential A a = A 



^ST> = e^e^—, = v?C-- - "-^) (18) 



dA x a ,dAP dA a 



dx 5 y dx a dxP 
Substitution of JTJJ and(THJ into (JT^J yields 



d P a , ,dA° dA a , mr dAi dA a ^ 



where (3 
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If we use in <|19|) the antisymmetric tensor F lk (e.g., [2]) 

nk dA k dA 

F = -dx~ ~ 8x7 (20) 

the relation between contra- (A a ) and co- variant (A a ) vectors (A a = — A a ), 
introduce the space-time interval 



ds = dty/l - f3 a /3 a = dty/1 - (3 2 

and the four-velocity 
we get 

^f- = qF ak u k = ~qF ka u k (21) 
ds 

As a next step, we find the zeroth components of Ecis. H21|) . Using the special- 
relativistic identity for the momentum pi = miii 

PiP 1 = m 2 
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wc find 

dp dp a a dp a 

Po— = -p a —=p a — (22) 
ds ds ds 

Upon insertion of (|22|l into l|21|l we obtain 

p a -^ = qp a [F^u (3 + F a0 u } (23) 

On the other hand, since F lk = -F ki { F 00 = F QQ = 0) 

p a u F al3 = 
Hence from J23 and (jSHJ) follows that 
drP 

_t_ = qu ^ Fa a = qF o aUa = qF 0i u . ( 24 ) 
ds 

Adding Eqs.J23J an( l I2D an d using the definition of F lk , Ea. l(2l7|) . we arrive at 
the equation of motion in the contra-variant form: 

The respective co- variant form follows from raising and lowering indices in (|25l) : 
d Pi r, k , dA h dAi k 



Now we reduce these equations to a form similar to the condition defining an 
irrotational flow in fluid mechanics. To this end we rewrite (|25[) and (|26Jl in the 
following form 

u k [^-(mu l + qA 1 ) - ^-{qA k )\ = 
OXk OXi 

A^imu, + qA) - ^-(qA k )] = (27) 
and add to the third term the identity 



Uk du k _ uk duk = 1 d , UkU k\ _ 
dxi dx l 2 dxi 



As a result, we get 



Pi Pi 

u k [—(mu l + qA 1 ) - g^(mu k + qA k )] = (28) 

or equivalently 

u k [^-r(mu i + qA i )-—(mu k + qA k )] = (29) 
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The expressions in square brackets represent a four-curl of the four- vector mui + 
qAi (or mu 1 + qA l ). Both equations are identically satisfied if this four-curl is 0. 
Once again, this can be interpreted as the fact that the respective vector field 
is irrotational, that is the four- vector mu + qA (here we use notation a for a 
four- vector) is the four-gradient of a scalar function, say — S 

mu 1 + qA l = -— (30) 

OXi 

dS 

mui + qAi = (31) 



This scalar function S ( a "potential function") is the classical relativistic ac- 
tion, and our choice of the sign is dictated by the consideration that expressions 
l|3U|) must become the expressions for the momentum and energy in the non- 
relativistic limit. 



To find the explicit expression for S we integrate Eq. (|30|l [ or l|31|) ] and obtain: 
S=- (mu 4 + qA^dxi = - / (m + A l Ui )ds (32) 



where a and b are points on the world line of the particle, ds — (dx l dxi) 1 ^ 2 , and 
Ui = dxi/ds. Expression i|32[) coincides ( as it should be) with the conventional 
definition of the action (introduced on the basis of considerations not connected 
to the second law of Newton). It is interesting to note that in a conventional 
approach to the action, the term A % dxi "cannot be fixed on the basis of general 
considerations alone" Here however this term is "fixed" by the very nature 
of the equations of motion. 

Eqs. H3()|l and 1)31(1 yield the determining PDE for the function S (the relativistic 
Hamilton-Jacobi equation for a charged particle in the electromagnetic field) if 
we eliminate Ui and u % from this equations with the help of the identity UiU 1 = 1: 

( S + qAl){ ^ + qA) = * = °' 1; 2 ' 3 (33) 

where we have to retain ( in the classical region) only one sign, either plus or 
minus. 1 



x We would like to point out that a unified way to describe wave and particle phenomena 
inherent to the Hamilton-Jacobi equation ( which was the main motivation of Hamilton) is 
conventionally demonstrated by comparing it and the eikonal equation and by showing that 
they are identical. On the other hand, as we show in the Appendix, there exists a simple 
way to do that directly from the Hamilton-Jacobi equation without resorting to the eikonal 
equation. 
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The usual way to derive the equations of motion (|^3f> [or from the action, 
Eq.J2U is to vary it. Here we follow the well-known procedure of reducing the 
integration of the partial differential equation of the first order to the integration 
of a system of the respective ordinary differential equations 1 . In particular, 
given the Hamilton- Jacobi equation ll-'i/it we derive (|25|) . To this end we subject 
action S to small perturbations 5S 

S = S +SS (34) 

and find the equation governing these perturbations. Here Sq must satisfy the 
original unperturbed Hamilton- Jacobi equation (|33[1 , and SS <C Sq. 

Upon substitution of (|34|) into (|33|) we get with accuracy to the first order in 
SS 

+ qA^SS + (|* + qAJ^-SS = (35) 
axi ox 1 ox 1 axi 

or equivalently 

pjQ Pi 

fe + ^)-*S = (36) 

Equation 1)36(1 is a quasi-linear first-order PDE whose characteristics are given 
by the following equations 

dx _ dx a 

dSo/dx° + qA dS /dx a + qA a { ' 

Here the repeated indices do not represent summation, and a = 1,2,3. It is 
immediately seen that the characteristics of linearized Hamilton-Jacobi equation 
l|37l) are the four- velocity u l : 

tt i = I(^ +9 A«) (38) 

' OXi 
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Inversely, these characteristics are the solutions of the equations of motion writ- 
ten in a form of the second law of Newton. To demonstrate that we divide both 
sides of l(38ll by ds, use Eqs. I(3Q|I . ((3I|I and the fact that d/ds ~ Ukd/dxk and 
obtain 

du l I ,dS Q . , d ,dS 
mc— = — + qA k — — ^ + qA % ) = 
as m Ox K Ox k OXi 

l^ +q A k )[— +q A l ) + q— -q— } = 

m dx k dxk dxi dxi dxi 

1 dS + qAk ^_d_^dSo + qAk ^ _ q( dA*_ _ dA* . _ 
m dx k dxi dxk dxi dxu 

1 d , ks 1 ,9A k dA\ - k 

^ 7 ^(u k u k ) + —qu k (- ) = qUk F tk (39) 

2m OXi m OXi ox k 
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that is the second law of Newton, Eq. 1251) 



Now we return to the linearized equation (|36fl which we rewrite in the identical 
form 

mu. l —SS=—SS = Q (40) 
oxi as 

This means that SS — const along a certain world line, singled out of a contin- 
uous set of possible world lines according to this condition. Without any loss of 
generality we can take the above const = 0. 

Thus on one hand, for a specific function S satisfying the Hamilton- Jacobi equa- 
tion the respective perturbations SS = const. On the other hand, according to 
Eq. (|32(l the action S is defined on a set of all possible world lines connecting 
world points a and b. This means that perturbations SS correspond to per- 
turbations of all these world lines. However, only for one of these world lines 
SS = const, according to l|40|) . To determine this constant we take into account 
that at the fixed world points a and b the world lines are also fixed, that is the 
respective perturbations SS = at these points. If we apply condition Ea. (|4*0")l 
to the action S, Eq. (|3f)|> . the former would "choose" out of all possible world 
lines the only one satisfying that condition, that is we arrive at the classical 
principle of least action. 

6 

(W + qA l )d Xl = (41) 



Finally, we demonstrate in an elementary fashion how the same technique of 
transforming the equations of motion in the Newtonian form to the Hamilton- 
Jacobi equation can be applied to a motion of a charged particle in general 
relativity. The equations of motion of a charged particle in gravitational and 
electromagnetic field are 3 . 

M(u 1 -^ + r kl u k u l ) = qg m F mk u k (42) 

where 

T i = 1 im ( d 9rnk ggmi _ dg ki 

kl 2 9 [ dx l dx k dx m > 



is the Ricci tensor. The expression T l kl u k u l is significantly simplified according 
to the following identities: 



r kl uu=u-[u(—^- 9mk -^ r )-u 9lm -^-ug 



l r I d 9 im , i k ( dg im dg u 1 [9 , Og™ , k dg kl 
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u dg lm u m 
1 dx 1 


-9 m 


du m 
dx 1 ' 


-5- 




,dg k[ u k 

( g xl 9M 


du k , 
da^' 




dx 1 



l.k d 9kl _ rdUi ii k du k 



dxi dx 1 ® dx 1 

j ik du k 1 9u\ 

dx 1 dx 1 2 dxi dxi 

- u l ^l + u l g tk ^± - u l — = -u'— + u l g ik (^^- - dUl ) 
9a;' 9x ; ctai da;' 9a;' da;' 

If we substitute this result into l|42|l and use the expression <|20|) for F , we 
obtain 

g^u^—^Muk + qA k ) - g^(M Ul + g^)] = (43) 
Equation (|4^|l is identically satisfied if we set 



Mu k + qA k = ( 44 ) 



where S is the action and we use the negative sign, representing a conventional 
choice of positive energies in classical mechanics. Raising and lowering the 
indices in l|44(l . expressing the respective 4- velocities u k and u k in terms of 
dS/dx k , and using the identity g lk UiU k = 1, we arrive at the Hamilton- Jacobi 
equation: 



We have shown that with the help of elementary operations one can arrive at 
the Hamilton- Jacobi equation from the phenomenological second law of Newton, 
without using of a priori defined action 5*. The latter arises in a natural way 
as a consequence of the existence of the "irrotational" solutions to the second 
law of Newton. The procedure follows from the fact that for forces determined 
by the potential energy U, the second law of Newton has a symmetry which 
allows us to reduce it to an antisymmetric form analogous to the form observed 
in the potential flow of an ideal fluid. This form lends itself to the introduction 
of a certain potential function (action), whose gradient represent a generalized 
momentum. Upon introduction of the action back into the second law of Newton 
we arrive at the Hamilton- Jacobi equation. We also show in an elementary 
fashion the wave-particle duality inherent in the Hamilton-Jacobi equation. As 
a result of this derivation we obtain the energy- momentum operators of quantum 
mechanics. 
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1 Appendix 



Let us consider a motion of a free relativistic particle of a mass m, whose 
Hamilton-Jacobi equation (|30[1 

dS OS , 
— -— — = m 

OX 1 OXi 

has a solution: 

S = - Pl x l (46) 

corresponding to the energy-momentum relation pip 1 — m 2 of the special rela- 
tivity. 

On the other hand, if we introduce the function 2 

S = hLnV (47) 

( where h is some constant having the dimension of the action S) the Hamilton- 
Jacobi equation l|46|) yields: 

d^dx~r mH (48) 

This equation admits the wave solution 

* = e~ lk ^ 3 (49) 
where k^iuj, k) is the wave four- vector. The respective dispersion relation is 

kjW = m 2 (50) 

Thus, on one hand the Hamilton-Jacobi equation describes a free particle with 
the momentum- energy four vector p l (E,p) and, on the other hand, the same 
equation describes a monochromatic wave with the wave four-vector k 1 

Returning back to function S ( Ea l48fl we obtain from (|49|l 

S = —hikjX 3 (51) 

Comparing Eas. (|46|l and (|51() we obtain the well-known relation between the 
wave four-vector and the four-momentum vector: 

*> - 1 < 52) 

Moreover, relation (|52() allows us to reinterpret the four-momentum vector in 
terms of the eigenvalues of the certain differential operators. In fact, from l|49|l 
and (|S2I> follows 

Mg^j = Pi* (53) 

2 which was done for the first time by E.Schroedinger in his historical paper [1] on non- 
rclativistic quantum mechanics 
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which means that we arrive at the relations describing quantum-mechanical 
operators of momentum-energy. 
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